Rectangular function
ββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββ
top
The rectangular function (also known as the rectangle function, rect function, Pi function, Heaviside Pi function,cite-ref-1[1] gate function, unit pulse, or the normalized boxcar function) is defined ascite-ref-wolfram-2-0[2]
rect β‘ ( t a ) = Ξ ( t a ) = { 0 , if | t | > a 2 1 2 , if | t | = a 2 1 , if | t | < a 2 . {\displaystyle \operatorname {rect} \left({\frac {t}{a}}\right)=\Pi \left({\frac {t}{a}}\right)=\left\{{\begin{array}{rl}0,&{\text{if }}|t|>{\frac {a}{2}}\\{\frac {1}{2}},&{\text{if }}|t|={\frac {a}{2}}\\1,&{\text{if }}|t|<{\frac {a}{2}}.\end{array}}\right.}
Its periodic version is called a rectangular wave.
Contents
β’ History
β’ See also
β’ References
ββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββββ
History
The rect function has been introduced 1953 by Woodwardcite-ref-6[6] in "Probability and Information Theory, with Applications to Radar"cite-ref-7[7] as an ideal cutout operator, together with the sinc functioncite-ref-8[8]cite-ref-9[9] as an ideal interpolation operator, and their counter operations which are sampling (comb operator) and replicating (rep operator), respectively.
Relation to the boxcar function
The rectangular function is a special case of the more general boxcar function:
rect β‘ ( t β X Y ) = H ( t β ( X β Y / 2 ) ) β H ( t β ( X + Y / 2 ) ) = H ( t β X + Y / 2 ) β H ( t β X β Y / 2 ) {\displaystyle \operatorname {rect} \left({\frac {t-X}{Y}}\right)=H(t-(X-Y/2))-H(t-(X+Y/2))=H(t-X+Y/2)-H(t-X-Y/2)}
where H ( x ) {\displaystyle H(x)} is the Heaviside step function; the function is centered at X {\displaystyle X} and has duration Y {\displaystyle Y} , from X β Y / 2 {\displaystyle X-Y/2} to X + Y / 2. {\displaystyle X+Y/2.}
Fourier transform of the rectangular function
The unitary Fourier transforms of the rectangular function arecite-ref-wolfram-2-1[2] β« β β β rect β‘ ( t ) β
e β i 2 Ο f t d t = sin β‘ ( Ο f ) Ο f = sinc β‘ ( Ο f ) = sinc Ο β‘ ( f ) , {\displaystyle \int _{-\infty }^{\infty }\operatorname {rect} (t)\cdot e^{-i2\pi ft}\,dt={\frac {\sin(\pi f)}{\pi f}}=\operatorname {sinc} (\pi f)=\operatorname {sinc} _{\pi }(f),} using ordinary frequency f, where sinc Ο {\displaystyle \operatorname {sinc} _{\pi }} is the normalized formcite-ref-10[10] of the sinc function and 1 2 Ο β« β β β rect β‘ ( t ) β
e β i Ο t d t = 1 2 Ο β
sin β‘ ( Ο / 2 ) Ο / 2 = 1 2 Ο β
sinc β‘ ( Ο / 2 ) , {\displaystyle {\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }\operatorname {rect} (t)\cdot e^{-i\omega t}\,dt={\frac {1}{\sqrt {2\pi }}}\cdot {\frac {\sin \left(\omega /2\right)}{\omega /2}}={\frac {1}{\sqrt {2\pi }}}\cdot \operatorname {sinc} \left(\omega /2\right),} using angular frequency Ο {\displaystyle \omega } , where sinc {\displaystyle \operatorname {sinc} } is the unnormalized form of the sinc function.
For rect β‘ ( x / a ) {\displaystyle \operatorname {rect} (x/a)} , its Fourier transform is β« β β β rect β‘ ( t a ) β
e β i 2 Ο f t d t = a sin β‘ ( Ο a f ) Ο a f = a sinc Ο β‘ ( a f ) . {\displaystyle \int _{-\infty }^{\infty }\operatorname {rect} \left({\frac {t}{a}}\right)\cdot e^{-i2\pi ft}\,dt=a{\frac {\sin(\pi af)}{\pi af}}=a\ \operatorname {sinc} _{\pi }{(af)}.}
Relation to the triangular function
We can define the triangular function as the convolution of two rectangular functions:
t r i ( t / T ) = r e c t ( 2 t / T ) β r e c t ( 2 t / T ) . {\displaystyle \operatorname {tri(t/T)} =\operatorname {rect(2t/T)} *\operatorname {rect(2t/T)} .\,}
Use in probability
Viewing the rectangular function as a probability density function, it is a special case of the continuous uniform distribution with a = β 1 / 2 , b = 1 / 2. {\displaystyle a=-1/2,b=1/2.} The characteristic function is
Ο ( k ) = sin β‘ ( k / 2 ) k / 2 , {\displaystyle \varphi (k)={\frac {\sin(k/2)}{k/2}},}
and its moment-generating function is
M ( k ) = sinh β‘ ( k / 2 ) k / 2 , {\displaystyle M(k)={\frac {\sinh(k/2)}{k/2}},}
where sinh β‘ ( t ) {\displaystyle \sinh(t)} is the hyperbolic sine function.
Rational approximation
The pulse function may also be expressed as a limit of a rational function:
Ξ ( t ) = lim n β β , n β ( Z ) 1 ( 2 t ) 2 n + 1 . {\displaystyle \Pi (t)=\lim _{n\rightarrow \infty ,n\in \mathbb {(} Z)}{\frac {1}{(2t)^{2n}+1}}.}
Demonstration of validity
First, we consider the case where | t | < 1 2 . {\textstyle |t|<{\frac {1}{2}}.} Notice that the term ( 2 t ) 2 n {\textstyle (2t)^{2n}} is always positive for integer n . {\displaystyle n.} However, 2 t < 1 {\displaystyle 2t<1} and hence ( 2 t ) 2 n {\textstyle (2t)^{2n}} approaches zero for large n . {\displaystyle n.}
It follows that: lim n β β , n β ( Z ) 1 ( 2 t ) 2 n + 1 = 1 0 + 1 = 1 , | t | < 1 2 . {\displaystyle \lim _{n\rightarrow \infty ,n\in \mathbb {(} Z)}{\frac {1}{(2t)^{2n}+1}}={\frac {1}{0+1}}=1,|t|<{\tfrac {1}{2}}.}
Second, we consider the case where | t | > 1 2 . {\textstyle |t|>{\frac {1}{2}}.} Notice that the term ( 2 t ) 2 n {\textstyle (2t)^{2n}} is always positive for integer n . {\displaystyle n.} However, 2 t > 1 {\displaystyle 2t>1} and hence ( 2 t ) 2 n {\textstyle (2t)^{2n}} grows very large for large n . {\displaystyle n.}
It follows that: lim n β β , n β ( Z ) 1 ( 2 t ) 2 n + 1 = 1 + β + 1 = 0 , | t | > 1 2 . {\displaystyle \lim _{n\rightarrow \infty ,n\in \mathbb {(} Z)}{\frac {1}{(2t)^{2n}+1}}={\frac {1}{+\infty +1}}=0,|t|>{\tfrac {1}{2}}.}
Third, we consider the case where | t | = 1 2 . {\textstyle |t|={\frac {1}{2}}.} We may simply substitute in our equation:
lim n β β , n β ( Z ) 1 ( 2 t ) 2 n + 1 = lim n β β , n β ( Z ) 1 1 2 n + 1 = 1 1 + 1 = 1 2 . {\displaystyle \lim _{n\rightarrow \infty ,n\in \mathbb {(} Z)}{\frac {1}{(2t)^{2n}+1}}=\lim _{n\rightarrow \infty ,n\in \mathbb {(} Z)}{\frac {1}{1^{2n}+1}}={\frac {1}{1+1}}={\tfrac {1}{2}}.}
We see that it satisfies the definition of the pulse function. Therefore,
rect β‘ ( t ) = Ξ ( t ) = lim n β β , n β ( Z ) 1 ( 2 t ) 2 n + 1 = { 0 if | t | > 1 2 1 2 if | t | = 1 2 1 if | t | < 1 2 . {\displaystyle \operatorname {rect} (t)=\Pi (t)=\lim _{n\rightarrow \infty ,n\in \mathbb {(} Z)}{\frac {1}{(2t)^{2n}+1}}={\begin{cases}0&{\mbox{if }}|t|>{\frac {1}{2}}\\{\frac {1}{2}}&{\mbox{if }}|t|={\frac {1}{2}}\\1&{\mbox{if }}|t|<{\frac {1}{2}}.\\\end{cases}}}
Dirac delta function
The rectangle function can be used to represent the Dirac delta function Ξ΄ ( x ) {\displaystyle \delta (x)} .cite-ref-0-11-0[11] Specifically, Ξ΄ ( x ) = lim a β 0 1 a rect β‘ ( x a ) . {\displaystyle \delta (x)=\lim _{a\to 0}{\frac {1}{a}}\operatorname {rect} \left({\frac {x}{a}}\right).} For a function g ( x ) {\displaystyle g(x)} , its average over the width a {\displaystyle a} around 0 in the function domain is calculated as,
g a v g ( 0 ) = 1 a β« β β β d x g ( x ) rect β‘ ( x a ) . {\displaystyle g_{avg}(0)={\frac {1}{a}}\int \limits _{-\infty }^{\infty }dx\ g(x)\operatorname {rect} \left({\frac {x}{a}}\right).} To obtain g ( 0 ) {\displaystyle g(0)} , the following limit is applied,
g ( 0 ) = lim a β 0 1 a β« β β β d x g ( x ) rect β‘ ( x a ) {\displaystyle g(0)=\lim _{a\to 0}{\frac {1}{a}}\int \limits _{-\infty }^{\infty }dx\ g(x)\operatorname {rect} \left({\frac {x}{a}}\right)} and this can be written in terms of the Dirac delta function as, g ( 0 ) = β« β β β d x g ( x ) Ξ΄ ( x ) . {\displaystyle g(0)=\int \limits _{-\infty }^{\infty }dx\ g(x)\delta (x).} The Fourier transform of the Dirac delta function Ξ΄ ( t ) {\displaystyle \delta (t)} is
Ξ΄ ( f ) = β« β β β Ξ΄ ( t ) β
e β i 2 Ο f t d t = lim a β 0 1 a β« β β β rect β‘ ( t a ) β
e β i 2 Ο f t d t = lim a β 0 sinc β‘ ( a f ) . {\displaystyle \delta (f)=\int _{-\infty }^{\infty }\delta (t)\cdot e^{-i2\pi ft}\,dt=\lim _{a\to 0}{\frac {1}{a}}\int _{-\infty }^{\infty }\operatorname {rect} \left({\frac {t}{a}}\right)\cdot e^{-i2\pi ft}\,dt=\lim _{a\to 0}\operatorname {sinc} {(af)}.} where the sinc function here is the normalized sinc function. Because the first zero of the sinc function is at f = 1 / a {\displaystyle f=1/a} and a {\displaystyle a} goes to infinity, the Fourier transform of Ξ΄ ( t ) {\displaystyle \delta (t)} is
Ξ΄ ( f ) = 1 , {\displaystyle \delta (f)=1,} means that the frequency spectrum of the Dirac delta function is infinitely broad. As a pulse is shorten in time, it is larger in spectrum.
See also
References
cite-note-11. β citerefwolfram-research2008Wolfram Research (2008). "HeavisidePi, Wolfram Language function". Retrieved October 11, 2022.
cite-note-wolfram-22. β reference-mathworld-rectangle-functionciterefweissteinWeisstein, Eric W. "Rectangle Function". MathWorld.
cite-note-77. β citerefwoodward1953Woodward, Philipp M (1953). Probability and Information Theory, with Applications to Radar. Pergamon Press. p. 29.
cite-note-1010. β Wolfram MathWorld, https://mathworld.wolfram.com/SincFunction.html
cite-note-0-1111. β citerefkharebutolarajora2023Khare, Kedar; Butola, Mansi; Rajora, Sunaina (2023). "Chapter 2.4 Sampling by Averaging, Distributions and Delta Function". Fourier Optics and Computational Imaging (2nd ed.). Springer. pp. 15β16. doi:10.1007/978-3-031-18353-9. ISBN 978-3-031-18353-9.